025921 -1 9-T 


CNAIA-CR-18B057) [TASK REPORTS ON 
DEVELOPING TECHNIQUES FOR SCATTERING PY 3 
CwSoSITE STRUCTURES ANT TO GENERATE NcW 
SOLOTIuNS in diffraction theory using high, 

OP.OSR BOUNDARY CONDITIONS! Semiannual 


R 


G3/ 2 7 


N9 1-19305 
— THRU — 
N91-I9307 
Unci as 
0002651 


SEMI-ANNUAL REPORT FOR NASA GRANT NAG-2-541 


Author: John L. Volakis 


February 1991 


National Aeronautics and 
Space Administration 
Ames Research Center 
Moffett Field, CA 94035 
Grant NAG-2-541 



THE UNIVERSITY OF MICHIGAN 

Radiation Laboratory 
Department of Electrical Engineering 
and Computer Science 
Ann Arbor, Michigan 48109-2122 
USA 


025921 -1 9-T 


SEMI-ANNUAL REPORT 
FOR 

NASA Grant NAG-2-541 
NASA Technical Monitor: Alex Woo 


Institution : The Radiation Laboratory 

Department of Electrical Engineering 
and Computer Science 
The University of Michigan 
Ann Arbor, Ml 48109-2122 


Period Covered : 


September 1990 to February 1991 


Report Title : 


Semi-Annual Report 


Principal Investigator : 


John L. Volakis 
Telephone: (313) 764-0500 


Preface 


This semi-annual report describes our progress during the 
period from September 1990 to February 1991. Two specific tasks 
are described and each should be read independently. That is, figure 
and reference numbering is consecutive only within the description 
of the task. As can be expected the progress reports are very brief 
and the reader should refer to the referenced technical reports for 
detailed coverage. A total of 23 journal articles, 17 conference 
publications and 18 technical reports have been written since the 
beginning of the Grant. A complete list of these is given in pp. 40- 
45. 
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ABSTRACT 


Validated results are presented for the new 3D body of revolution finite element- 
boundary integral code. As usual, a Fourier Series expansion of the vector electric and 
magnetic fields is employed to reduce the dimensionality of the system and the exact 
boundary condition is employed to terminate the finite element mesh. The mesh 
termination boundary is chosen such that it leads to convolutional boundary operators for 
low 0(n) memory demand. Improvements of this code are discussed along with the 
proposed formulation for a full 3D implementation of the finite element-boundary integral 
method in conjunction with a CGFFT solution. 

OBJECTIVE 

The objective of this task is to develop innovative techniques and related software 
for scattering by three dimensional composite structures. The proposed analysis is a hybrid 
finite element-boundary integral method formulated to have an O(n) memory demand. This 
low storage is achieved by employing the FFT to evaluate all boundary integrals and by 
resorting to an iterative solution algorithm. Particular emphasis in this task is the 
generation of software applicable to airborne vehicles and the validation of these by 
comparison with measured and other reference data. Because the approach is new, a step 
by step development procedure has been proposed over a three-year period. During the 
first year the technique was developed and implemented for two-dimensional composite 
structures. Support software for the two-dimensional analysis such as pre- and post- 
processor routines were developed during the second year and a formulation was also 
developed and implemented for three-dimensional bodies of revolution. Finally, during the 
third year, we will develop, implement, and test the method for arbitrary three dimensional 
structures. 
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BACKGROUND 

Interest in three-dimensional (3-D) methods has increased in recent years, however, 
the associated demands in computation time and storage are often prohibitive for electrically 
large 3-D bodies. Vector and concurrent (i.e. hypercube, connection, etc.) computers are 
beginning to alleviate the first of these demands, but a minimization of the storage 
requirements is essential for treating large structures. 

The traditional Conjugate Gradient Fast Fourier Transform (CGFFT) method [1 ] - 
[4] is one such frequency domain solution approach which requires O(n) storage for the 
solution on n equations. This method involves the use of FFTs whose dimension equals 
that of the structure under consideration [5] - [7] and, therefore, demands excessive 
computation time when used in an iterative algorithm. Also, the standard CGFFT requires 
uniform rectangular gridding that unnecessarily includes the impenetrable portions of the 
scatterer. With these issues in mind, a new solution approach is proposed for solving 
scattering problems. The proposed method will be referred to as the Finite Element- 
Conjugate Gradient Fast Fourier Transform (FE-CGFFT) method. 

During last year's effort the FE-CGFFT method was developed for two- 
dimensional scatterers where the finite element mesh was terminated at a rectangular box. 
Inside the box boundaries, Helmholtz equation is solved via the finite element method and 
the boundary constraint is obtained by an appropriate integral equation which implicitly 
satisfies the radiation condition. Along the parallel sides of the box, this integral becomes a 
convolution and is, therefore, amenable to evaluation via the FFT. The dimension of the 
required FFT in this hybrid method is one less than the dimensionality of the structure thus, 
making it attractive for 3-D simulations. Also, because it incorporates the finite element 
method, the FE-CGFFT formulation remains valid regardless of the structure's geometry 
and material composition. 

The proposed method described in the University of Michigan Report 02592 1-6-T 
(see also [8]) is similar to the moment method version developed by Jin [9]. Jin's method 
was in turn based on work published in the early 70’s by McDonald and Wexler [10] who 
introduced an approach to solve unbounded field problems. The proposed method is also 
similar to other methods (a few of which will be mentioned here), neither of which 
provides a storage reduction comparable to the proposed FE-CGFFT method. The 
unimoment method [11] uses finite elements inside a fictitious circular boundary and an 
eigenfunction expansion to represent the field in the external region. The coefficients of the 
expansion are then determined by enforcing field continuity at the finite element (FE) mesh 
boundary. The coupled finite element-boundary element method [12] uses the finite 
element method within the boundary and the boundary element method to provide the 
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additional constraint at the termination of the mesh. Unlike the proposed method, the 
solution in [12] was accomplished by direct matrix inversion (as in [9]), and the outer mesh 
boundary is not rectangular to take advantage of the FFT for the evaluation of the boundary 
integrals. 


P-RQGRESS 

Part of our efforts in this task were devoted to debugging and validating the three 
dimensional body of revolution (BOR) code developed in the previous months. The 
analysis associated with this code is described in the U of M technical report 02592 1-18-T 
where we also include validation data obtained over the past two months. Some of these 
are shown in figures 1-3 and refer to an ogive, a circular cylinder and a sphere. 
Unfortunately, it was found that as the bodies became larger the system's condition 
deteriorates and this was traced to the pulse basis formulation employed for the 
discretization of the boundary. Through several tests we have now shown that A 
Galerkin's linear basis formulation will correct the convergence difficulties. For example, 
this formulation was already employed in solving large systems (with more than 120,000 
unknowns) associated with the scattering by frequency selective surfaces (FSS) and large 
Mates. As shown in figure 4-5, the Galerkin's formulation with linear basis permitted a 
solution of this system in less than 70 iterations! In comparison, the pulse basis-point 
matching formulation would require several thousand iterations before reaching 
convergence. Consequently, we are in the process of incorporating the Galerkin's linear 
basis formulation into our existing 2D and 3D BOR codes. Further, it was found that the 
ech area converges much sooner than the mean square error and permitted us to speed-up 
solution time. 

During this last quarter we also began the development of the proposed finite 
element formulation for general non-symmetric inhomogeneous bodies. The basic discrete 
elements in this case are tetrahedra in conjunction with edge-based expansion functions. 
The associated finite element formulations is described in Appendix A and we are now in 
the process of implementing it. Initially, the finite element mesh will be terminated by a 
fictitious absorbing layer whose dielectric parameters were determined by a minimization of 
the reflection coefficient over the entire range of incidence angles. For a three layer 
coating, each of thickness 0.05 wavelengths, it was found that their respective dielectric 
properties to minimize the reflection coefficient over all angles of incidence are 
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e r i =(-0.1249205, -1.731605), |j. r i=(- 1.03 1792; -0.1039932) 
e r 2=(0.040699530, 0.1750280), p r i=(0.3155941, 0.3190330) 
e r 3 =(- 1.278644, 0.9625375), |X r 3=(-0. 17213 15, -5.389832) 

The corresponding plot of the reflection coefficient as a function of incidence is given in 
figure 7 along with scattering patterns based on the proposed termination model. As seen, 
for the chosen fictitious absorber the reflection coefficient is less than one percent for 0 up 
to 62 degrees and less than 2 percent for 0 up to 77 degrees. For the same error criteria, 
the corresponding angles associated with the second order Pade ABC are 35 and 41 
degrees, respectively. The fictitious ABC has, therefore, a substantially better performance 
over the existing ABCs, and its effectiveness will be examined further in the next few 
months. 

The three dimensional finite element meshes required in the analysis will be 
generated by SDRC IDEAS and we have already began to develop the software for 
transforming the output of this commercial package to the input files of our analysis codes. 
Similar drivers were already developed for the two dimensional code which was developed 
last year. 

Finally, during this period we performed extensive testing of the two-dimensional 
code and have in the progress developed several new pre- and post- processing algorithms 
for this code. Two of the new geometries (see fig. 8 and 11) whose scattering was 
computed with our 2D finite element - CGFFT code are displayed in figures 9, 1 1 and 1 2. 
These represent airfoil configurations, one of which is coated with a dielectric material. 

CONCLUSIONS 

The project continues to evolve in accordance with our original plan and schedule. 
Most importantly, so far, our expectation of the finite element CGFFT formulation have 
been realized and we are, therefore, pleased with its performance for the intended 
applications. 

TRANSITIONS 

All of our efforts in the next six months will be devoted to the development of the 
3D finite element boundary integral code for arbitrary structures. In the immediate future 
we will also pursue improvements for our existing codes primarily directed at speeding the 
convergence of the CG or BiCG algorithm. 
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Body of Revolution enclosed in a Fictitious finite length cylinder. Area between 
the Scatterer and the fictitious cylinder is discretized via the finite element method. 
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Figure 1: TM and TE bistatic scattering 
pattern from a perfectly conducting sphere of 
radius for axial incidence. 
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Figure 4: Scattering by a 12 x 12 FSS array; comparison of the exact solution 
(solid line) with an approximate result obtained by truncating the infinite FSS. 
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FSS ARRAY SCATTERING 
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Figure 5: Scattering by a 24 x 24 array; comparison of the exact solution 
(solid line) with an approximate result obtained by truncating the infinite FSS. 
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Figure 7: Evaluation of the new fictitious ABC (a) geometry (b) reflection coefficient 

plot (c) H z field on a PEC Circular Cylinder (d) E z field at a distance X/10 from the 
surface of a PEC Circular Cylinder. 
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case 2 : Coated Trailing Edge 

±±V l-(x/2.5) 2 -2.5<x<0 

±0.8232 A(x) 0<x<2.5 

A(x) = V(l-(x/2.7182f) - 0.3926 

y out = ±0.81 16 B(x) 0.1<x<3.0 

y 0 ut=yin elsewhere 

B(x) = V l-[(x-.l)/3.1416] 2 -.3846 

e r =2-j 1 between y in and y out 
FEM Mesh 
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Case 2: Coated trailing edge, (£=2-jl), H-pol 



Angle [deg] 



Angle [deg] 


Figure 9: E and H polarization scattering patterns for the configuration shown in 
o. Comparison of results between the FE-CGFFT and the moment method. 


figure 


16 




Perfectly conducting airfoil 

All dimensions in wavelengths. The airfoil section is made by 5 arcs: 



• OA : straight line 

• AB : circle of radius R2 = 7Aq and of center 02 

• BC : polynomial parametric equation 

• CD : polynomial parametric equation 

• DO : circle of radius Rl = 9Ao and of center 01 
The polynomial equation are given by: 

*(«) = * « (# ^ »(«) * Ej.,4, * 0 < a < 1 

* “d ^ for BC are «, *nd 4, for CD arc 

* 4.61149 = 1.53278 = >-1.62131 4i = -0.12563 

= -12.11403 4, = -3.22680 «i = 4.54389 5,= 0.30612 

«,* 8.88606 6, = 0.90615 s, = -3.32901 6,= 0.16113 

* = 0.44275 4«a 1.24623 o< a 0.35663 4 4 = -0.00343 

«• * —0.51440 6| * 0.46927 o« = -1.40984 b% a -0.51216 

0 *a -6.86720 4,= -0.41007 o« = -5.40756 4, a -0.23609 

Figure 10: Geometry of a PEC Airfoil whose scattering is given in figures 11 & 12. 
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Appendix A 


Finite element formulation for tetrahedral 
elements and edge-based expansion basis 

1 Derivation of finite element equations 

Let us consider a three dimensional inhomogeneous body occupying the 
volume V. In order to discretize the electric field E inside the body, we 
subdivide the volume V into a number of small tetrahedra, each occupying 
volume V e (e = 1,2, • ■ • , M ) with M being the total number of tetrahedral 
elements. Within each tetrahedron, the electric field satisfies the vector 
wave equation 

V x —V x E - k\tr E = 0 (1) 

/z r 

where /.i r is the permeability of the medium, e r is the medium permittivity 
and k a is the free space wave number. The next step is to expand the 
electric field within V e as 

E = ££*W‘ (2) 

where W® are edge-based vector basis functions and E e - denote the 
expansion coefficients of the basis, all defined within the volume V e . W| is 
tangential to the jth edge of the eth tetrahedron with zero tangential 
component along the other edges of the tetrahedral element. On 
substituting (2) into (1), we obtain 

E^'fvxivx Wj - ^e r w;) = 0 (3) 

In order to solve for the unknown expansion coefficients Ej, we take the dot 
product of (3) with W* and then integrate the resulting equation over the 
element volume V e (Galerkin’s technique). The wave equation thus reduces 
to 
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E E ) X W '-( Vx 7 Vx W i “ J dv = 0 (4) 

The first term in the integral of the above expression can be simplified by 
using basic vector identities. Since 

W'. Vx —V x W'j = V • —(V x W*) x Wf + —(V x W?) . (V x W') 

the divergence theorem can be readily applied to (4) resulting in the 
following expression: 


o = 


-(V x Wf) . (V x W*) - k\t r W, e .W' I dv 


jTl [ JV ‘ \Vr 


+ £ t (VxW ' )xW ') dS } (5) 

where S e denotes the surface enclosing V e . Using vector identities , (5) can 
be further simplified to yield the weak form of Maxwell’s equation: 

j^E' [ (—(V x W') . (V x W') - dv = jufi 0 <f W'.(n x H )ds 

j=l •'V* \f*T ) 


where n x H is the tangential magnetic field on the exterior dielectric 
surface. Equation (6) can be conveniently written in matrix form as 


Ml*:-] = [S'] 


where 


A h = / (-(VxW').(VxW;)-i:£ r W'W (S) 

JVt \ ft r J 

B' = jup 0 l W e .(n x H )ds (9) 

JSe 

On assembling all the M tetrahedral elements that make up the geometry, 
we obtain a system of equations whose solution yields the field components 
over the entire body. Therefore, summing over all M elements, we have 
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( 10 ) 


M M 

Emm = Eim 

e=l e=l 

which gives 

[A)[E] = [ B ] (11) 


where [A] is a N x N matrix with N being the total number of edges 
resulting from the subdivision of the body and [ E ] is a N x 1 column vector 
denoting the edge fields. Due to the continuity of the tangential component 
of the magnetic field at the interface between two dielectrics, an element 
face lying inside the body does not contribute to [ B] since the surface 
integrals over the faces of adjacent tetrahedra cancel each other. As a result, 
[B] is a column vector containing the tangential magnetic field only over 
the exterior surface of the body. Equation (11) can therefore be written as 

+ Aa,’E; = H, 

A l 3 E s + A it E t = 0 (12) 


where the subscript s denotes the edges on the surface and i represents the 
edges inside the body. It is thus readily seen that (11) relates the electric 
field inside and on the surface of the body to the on-surface tangential 
magnetic field. 


2 Basis functions 


Vector fields within tetrahedral domains in three dimensional space can be 
conveniently represented by expansion functions that are linear in the 
spatial variables and have either zero divergence or zero curl. The basis 
functions defined below are associated with the six edges of the tetrahedron 
and have zero divergence and constant curl. Assuming the four nodes and 
the six edges of a tetrahedron are numbered according to Table 1, the 
vector basis functions associated with the (7 — i)th edge of the tetrahedron 
are defined as 


W 


7—i — 


{ 7 + g 7 _, x r, r in the tetrahedron 
0, otherwise 


where i = 1, 2, ■ • • , 6 and f and g are constant vectors. On direct 
evaluation, it is readily seen that 


(13) 


22 



V-W, = 0 
VxW, = 2g, 


(14) 

(15) 


Since the complex scalar Ej in ( 2) is the projection of the electric field 
onto the jth edge of the tetrahedral element, 


w -e = 8 - 

1 r on jth edge 


(16) 


where <5jj is the Kronecker delta. Solving ( 13) and ( 16) for the unknown 
vectors yield[l] 


fz-, 


g7-i 


h-i 

6V r “ X r * 

^j^7— j6j 

6V 


(17) 

(18) 


where V is the volume of the tetrahedral element, = (r,- 2 — r fl )/6,- is the 
unit vector of the ith edge and 6, = |ri a — r,-, | is the length of the ith edge. 
All distances are measured with respect to the origin. 

Since there are two numbering systems, local and global, a unique global 
direction is defined (e.g., always pointing from the smaller node number to 
the larger node number) to ensure the continuity of n x E across all edges. 
This implies that ( 13) should be multiplied by (-1) if the local edge vector 
(as defined in Table 1) does not have the same direction as the global edge 
direction. Even though W, forces no conditions on the normal component 
of E, it has been shown[2] that the continuity of electric flux can be satisfied 
within the degree of approximation with the above formulation. Finally, 
since V • W, = 0 the electric field obtained through ( 2) exactly satisfies the 
divergence equation within the element, i.e. V ■ E = 0. Therefore, the finite 
element solution is free from contamination of spurious solutions[2]. 


3 Mesh termination 

Differential equation methods, such as finite elements, can only solve 
boundary value problems. Since electromagnetic problems are open 
boundary-infinite domain types, a means to truncate the solution domain 
to lie within a finite boundary must be found. On this boundary, a 
condition is enforced thus ensuring that the fields will obey the Sommerfeld 
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radiation condition at distances asymptotically far from the object. These 
absorbing boundary conditions (ABCs) have a significant advantage over 
the global methods of solving unbounded problems using finite elements in 
that they are local in nature. Due to this, the sparse matrix structure of 
the finite element formulation is retained. One disadvantage, however, is 
that ABCs are approximate and do not model the exterior field exactly. 

The objective of absorbing boundary conditions is to truncate the finite 
element mesh with boundary conditions that cause minimum reflections of 
an outgoing wave. These ABCs should provide small, acceptable errors 
while minimising the distance from the object of interest to the outer 
boundary. This minimal distance is required to reduce the number of 
unknowns in the problem for computational efficiency. A three dimensional 
vector boundary condition will be investigated here for terminating the 
finite element mesh of the body described in section 1.1. We begin with the 
Wilcox representation [3] of the electric field which has an expansion 


E(r) = 


e ~‘ k ' 


r 


n=0 7 


(19) 


From ( 19), we get 

V x E = | jkr x 1 + ^ £>1 1 E 


e-^ r ~ nA nt 

ir » 2 7»fl 

7 n=l 7 


( 20 ) 


where A n( = r x A n is the transverse component of A n and, for a vector F, 
DiF is given by 


DiF = 


1 

sinO 
( 1 

sinO 




8F T 

— - sindF* 

6 + 

k-^1 

86 


86 


* 


( 21 ) 


Using the recursion relation 
- 2jknA nt 


n(n - l)A n _i, t + D 4 A n _i 


where 



(D A e n + De A n ) 6 + (D At + D+ A n ) 4> 

}_ A e _ 2 cos0 8Aj 

86 sin 2 6 " sin 2 6 86 

J_?K _ _i_ , 2co£3K 

sin6 84> sin 2 0 n sin 2 6 8<j> 

and D is Beltrami’s operator [3], we can derive the representation correct to 
r~ 4 . Applying the recursion relation in ( 20) yields the desired relationship 
for the vector ABC: 


D 4 A n — 
A„ = 
D^ A n 


V x E = a(r)E + /?(r)Z) 4 E 


( 22 ) 


where 


o:(r) 

P(r) 



1 1 


2 jkr 2 (1 + 1/jkr) 


(23) 

(24) 


The ABC formulated above is applicable to spherical boundaries and hence 
would be storage intensive and numerically inefficient when used to 
terminate the mesh of long and thin geometries. It would be highly 
desirable to choose an outer boundary that conforms to the shape of the 
object. An approximate boundary condition based on the asymptotic 
representation of fields for a two dimensional scalar problem has already 
been derived[4j. It is the author’s intention to extend the derivation of the 
two dimensional scalar boundary condition to a three dimensional vector 
absorbing boundary condition for an arbitrary outer boundary. 


4 Solution of the finite element equations 

An inspection of ( 11) reveals that for an inhomogeneous body, there is no 
a priori information about the tangential magnetic field over the exterior 
surface of the body. Relation (11) therefore contains two unknown vectors, 
[E] and [B], and thus another condition is required involving the two 
variables to permit an evaluation of the fields inside and on the surface of 
the body. This condition relating the tangential electric field to the 
tangential magnetic field on the surface is provided by ( 22). Since the 
ABC in ( 22) refers to the scattered field, we can rewrite it as 


25 


VxEJ 

h; 


«(r)E; + fi(r)D 4 E’ 

— [a(r) E‘, + p(r)D 4 E’ 3 ] 

LOfl 

fc e : 


( 25 ) 


where /C = ^ [a(r) + /?(r)D 4 ] and the subscript s denotes the field on the 
surface and the superscript s represents the scattered field. Since the total 
field is a sum of the incident field and the scattered field, therefore from 
( 25), we obtain 


h; = /ce; 

h, -h;“ = k (e, - E; ne ) 


(26) 


Substituting ( 26) into ( 12) and simplifying gives 

(A„ - JC)E a + A ai Ei = H* nc -/CE* nc 

A.^E, + j4„E, = 0 (27) 


The above equation can thus be solved for the unknown electric fields both 
inside and on the surface of the body. 
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Tetrahedron Edge Definition 
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ABSTRACT 


Rigorous UTD (uniform Geometrical Theory of Diffraction) diffraction coefficients 
are presented for a coated convex cylinder simulated with generalized impedance boundary 
conditions. In particular, ray solutions are obtained which remain valid in the transition 
region and reduce uniformly to those in the deep lit and shadow regions. These involve 
new transition functions in place of the usual Fock-type integrals, characteristic to the 
impedance cylinder. A uniform asymptotic solution is also presented for observations in 
the close vicinity of the cylinder. As usual, the diffraction coefficients for the convex 
cylinder are obtained via a generalization of the corresponding ones for the circular 
cylinder. 


OBJECTIVE 

This task involves the use of higher order boundary conditions to generate new 
solutions in diffraction theory. In particular, diffraction coefficients will be developed for 
dielectric/magnetic layers and metal-dielectric junctions which are often encountered on 
airborne vehicles as terminations of coatings and conformal antennas. Solutions for both 
polarizations will be developed for fairly thick junctions and versatile computer codes will 
be written and tested. Creeping wave diffraction coefficients will be also developed for 
multilayered coated cylinders. 
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PROGRESS 


1 . Introduction 

The problem of scattering by a smooth convex impedance cylinder has received 
much attention. Wang [1,2] presented ray-optical solutions for the impedance and coated 
cylinders. His results are valid only in the deep lit and shadow regions and do not apply to 
the case where the observation point is in the transition region. Wait and Conda [3, 4] 
developed a solution which is valid in the transition region and for observation points on 
and off the surface. However, as pointed out by Pathak [5] it did not uniformly reduce to 
the ray solution [6, 7] exterior to the transition regions. Also, it is not valid on the portion 
of the surface in the transition region and these limitations were the primary motivation in 
Pathak’s work [5] for the perfectly conducting convex cylinder. Recently, Kim and Wang 
[8] presented a solution applicable to a coated cylinder that remained valid in the transition 
region. They employed a heuristic approach to obtain the numerical values of the resulting 
transition integral applicable to a coated cylinder. Their solution is uniform but is not 
applicable to the close vicinity of the cylinder. 

Here we develop a rigorous UTD solution of the diffraction by a coated cylinder 
simulated with generalized impedance boundary conditions. In addition, a uniform 
asymptotic solution is obtained which remains valid when the observation point is in close 
vicinity of the cylinder. An important aspect of the paper is also the use of second order 
generalized impedance boundary conditions (GIBC) for the simulation of the coating. 
Their derivation has already been given in [9] and [10] and are characterized by the 
inclusion of higher order field derivatives in their definition. Because of this they are less 
local which leads to an improved simulation (with respect to the standard impedance 
boundary condition - SIBC) of the coating in a manner analogous to the order of the 
highest derivative kept in the condition. Recently, they were successfully applied to a 
number of diffraction problems [11], [12] and have also been used in numerical 
simulations of multilayer coatings (see fig. 1) [13]. These applications provided a measure 
of the accuracy of the proposed GIBC and in particular accuracy criteria were derived in 
[13] for the second order conditions as a function of coating thickness and composition. 

The UTD solution to be presented here parallels that given by Pathak [5] for the 
circular perfectly conducting cylinder. However, in the case of the coated cylinder the 
resulting UTD expressions are in terms of Fock-type integrals whose efficient evaluation is 
of primary interest. In the following we first present the eigenfunction solution based on 
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the second order GIBC simulation of a circular coated cylinder. By employing Watson's 
transformation this is written in integral form which is then cast in a ray representation. 
They ray solution is subsequently generalized to the case of a general convex cylinder. 
Finally, the evaluation of the Fock-type integrals is discussed and some results are 
presented which validate the accuracy of the GIBC eigenfunction and ray solutions. In the 
process, we demonstrate the improved accuracy of the GIBC solution over the 
corresponding SIBC solution, and it is also shown how the presented UTD solution can be 
extended to treat multilayered coated cylinders. 

The details of the analysis are described in the report 025921 -13-T which was 
recently submitted to the sponsor. Below we only attach a few results which demonstrate 
the accuracy and utility of the derived formulae. 

2. Numerical Results 

The UTD expressions derived in the UM report 025921-17-T provide a complete 
set of equations for the computation of the total field in all regions of interest. Below, we 
present some calculated data which validate the accuracy of the derived expressions by 
comparison with data based on the moment method and eigenfunction solutions. 

In figure 2 the eigenfunction solutions based on the GIBC and SIBC simulations 
are compared with the exact for a coated cylinder and this clearly demonstrates the 
improved simulation (with respect to the standard impedance boundary conditions - SIBC) 
achieved with the second order GIBC. To show the validity of the UTD solution in the 
case of the convex cylinder, a special case of an elliptical cylinder is considered in figure 3. 
Data based on the moment method are compared with those obtained from the UTD 
solution in conjunction with the second order low and high contrast boundary conditions. 

Figure 4 verifies the asymptotic solution developed for the field point in the close 
vicinity of a convex cylinder. We remark, however, that the approximations used for the 
Hankel functions in the derivation of (42) and (43) become less accurate for some values of 
e r and |X r associated with lossless coatings, and this can be avoided by using more accurate 

approximations for the Hankel functions. Finally, figure 5 demonstrates the use of GIBC 
in simulating multilayer coatings by simply redefining the material constants a m and a m ' as 
discussed in [10, 13]. 

A difficulty in implementing the expressions derived in this paper was the 
evaluation of the Fock-type integrals G(x,q), gj (D) and g 2 (D) as well as determination of 
the zeros corresponding to (21). The Fock-type integrals were evaluated by employing the 
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method described in [16] and the zeros of (21) were determined using the routine given in 

[ 20 ]. 


Summary 

Rigorous ray solutions of the scattered fields were presented for a coated convex 
cylinder. These were developed in the context of the uniform geometrical theory of 
diffraction and specific expressions were given for the scattered fields in the lit, shadow 
and transition regions as well as for observations in the near vicinity of the cylinder. That 
is, UTD expressions were derived for all regions exterior to the coated cylinder. These are 
suited for engineering computations and are given in terms of the generalized Pekeris or 
Fock-type functions whose evaluation was efficiently performed via the Fourier 
Trapezoidal rule suggested by Pearson [16]. 

In comparison to the solution given by Kim and Wang [8], the ray representations 
given here are based on a second order generalized impedance boundary condition which 
permits the simulation of thin multilayered coating as demonstrated in the included 
examples. Also, in our implementation of the transition fields we employed a rigorous 
rather than a heuristic evaluation of the Fock-type integrals. Further, we have presented 
accurate field representations for observations on or near the vicinity of the coated cylinder 
and these can also be used for computing the radiated fields by a source or an aperture on 
the surface of the convex cylinder. 

REFERENCES 

[1] N. Wang, "Regge poles, natural frequencies, and surface wave resonance of a 
circular cylinder with a constant surface impedance," IEEE Trans. Antenna & 
Propagat., Vol. AP-30, Nov. 1982. 

[2] — , "Electromagnetic scattering from a dielectric coated circular cylinder," IEEE 
Trans. Antenna & Propagat., Vol. AP-33, Sept. 1985. 

[3] J.R. Wait and A.M. Conda, "Diffraction of electromagnetic waves by smooth 
obstacles for grazing angles," /. Res. Bur. Stand., Vol. 63D, pp. 181-197, 1959. 

[4] — , "Pattern of an antenna on a curved lossy surface," IRE Trans. Antennas & 
Propagat., Vol. AP-6, pp. 348-359, Oct. 1958. 


32 



[5] P.H. Pathak, "An asymptotic analysis of the scattering of plane waves by a smooth 
convex cylinder," Radio Sci., Vol. 14, pp. 419-435, 1979. 

[6] J.B. Keller, "Diffraction by a convex cylinder," IEEE Trans. Antenna & Propagat., 
Vol. AP-24, pp. 312-321, 1956. 

[7] — , "Geometrical theory of diffraction," J. Opt. Soc. Amer., Vol. 52, 
pp. 116-130, 1962. 

[8] H. Kim and N. Wang, "UTD solution for electromagnetic scattering by a circular 
cylinder with thin lossy coatings," IEEE Trans. Antenna & Propagat., Vol. AP-37, 
pp. 1463-1472, Nov. 1989. 

[9] T.B.A. Senior and J.L. Volakis, "Derivation and application of a class of 
generalized boundary conditions," IEEE Trans. Antenna & Propagat., Vol. AP-37, 
pp. 1566-1572, Dec. 1989. 

[10] M.A. Ricoy and J.L. Volakis, "Derivation of generalized transition/boundary 
conditions for planar multiple layer structures," Radio Sci., Vol. 25, pp. 391-405, 
July - August 1990. 

[11] — , "Diffraction by a symmetric material junction. Part I: general solution," 
submitted to IEEE Trans. Antenna & Propagat. 

[12] — , "Diffraction by a symmetric material junction. Part II: boundary conditions," 
submitted to IEEE Trans. Antenna & Propagat. 

[13] J.L. Volakis and H.H. Syed, "Application of higher order boundary conditions to 
scattering by multilayered coated cylinders," J. of Electromagnetic Waves and 
Applications, (in press) 

[14] J.J. Bowman, T.B.A. senior and P.L.E. Uslenghi, "Electromagnetic and Acoustic 
Scattering bv Simple Shapes ." Amsterdam: North-Holland, 1969. 

[15] N.A. Logan and K.S. Yee, " Electromagnetic Waves." edited by R.E.Langer, 
University of Wisconsin Press, Madison. 


33 


[16] L.W. Pearson, "A scheme for automatic computation of Fock-type integrals," IEEE 
Trans. Antenna & Propagat., Vol. AP-35, pp. 1111-1118, Oct. 1987. 

[17] O.A. Tuck, "A simple Pinion-trapezoid' rule,” Math. Comp., Vol. 21, 
pp. 239-241, 1967. 

[18] L. Felson and N. Marcuvitz, " Radiation and Scattering of Waves ." 

Prentice Hall, 1973. 

[19] R.G. Kouyoumjian and P.H. Pathak, "A uniform geometrical theory of diffraction 
for an edge in a perfectly conducting surface," Proc. IEEE, Vol. 62, 

pp. 1448-1461, 1974. 

[20] B.K. Singaraju, D.V. Giri and C.E. Baum, "Further development in the application 
of contour integration to the evaluation of the zero of analytical functions and 
relevent computer programs," Math. Note 24, March 1976. 


34 




35 


Total Field (dB) Total Field (dB) 



Observation Angle 4> (deg) 
(a) 
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Fig. 2: Bistatic H-polarization scattering pattern for a circular cylinder 

having b = 3A, (a) e r = 4, /i r = 1, 6 = 0.07A, p = 5X (b) e r = 8 
Hr = 1, 6 = 0.2A, /> = 3.05A. 
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Flg ‘ 3: Bistatic scattering pattern for an elliptical cylinder having a = 2A, 

b = 1A, p = 5A, 4>i = 0 (a) E-polarization, e r = 4, = 1, 

6 = 0.07A (b) H-polarization, e r = 8, (i r = 1, 6 = 0.2 A. 
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(a) 



Observation Angle $ (deg) 
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Pig ♦ 

Bistatic scattering pattern of a circular cylinder having b = 3A, 
p = 3.05A, <tn = 0 (a) E-polarization, c, = 4, p r = l, 6 = 0.07A 
(b) H-polarization, e, = 8, = 1, 6 = 0.2A. 
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Observation Angle <|> (deg) 

(b) 

Bistatic scattering pattern of a three- layer coated circular cylinder 
having b = 3A, = 0, e rl = 3-;'0.1, u* = 4-j0.3, ^3 = 7-; 1.5, 

firi = /irj — /ira = 1, = 0.01A, 63 = 0.02A, 63 = 0.03A (a) 

H-polarization, p = 5 A (b) E- polarization, p = 3.05A. 
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